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The Michelson Interferometer is a device used in many applications, but here 
it was used to measure small differences in distance, in the milli-inch range, 
specifically for defects in the Orbiter windows. In this paper, the method of 
using the Michelson Interferometer for measuring small distances is explained 
as well as the mathematics of the system. The coherence length of several 
light sources was calculated in order to see just how small a defect could be 
measured. Since white light is a very broadband source, its coherence length 
is very short and thus can be used to measure small defects in glass. After 
finding the front and back reflections from a very thin glass slide with ease 
and calculating the thickness of it very accurately, it was concluded that this 
system could find and measure small defects on the Orbiter windows. This 
report also discusses a failed attempt for another use of this technology as 
well as describes an area of promise for further analysis. The latter of these 
areas has applications for finding possible defects in Orbiter windows without 
moving parts. 
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Introduction 


The Michelson Interferometer is implemented in many different fields, 
like medicine, detection of gravitational waves, sensing, and many others. The 
application that we are concerned with involves using the interference proper- 
ties of light to find small differences in length. For this project, the Michelson 
Interferometer will be used to measure small defects in the Orbiter windows. 
As of now, there are instruments and tools that can detect defects to a few 
micro-inches, but it takes a large amount of time to do this. This project 
involves using the Michelson Interferometer so that scratches and defects on 
the surface of the Orbiter windows can be pinpointed and dealt with in a 
timely fashion. A defect in a window is of concern on the Orbiter windows 
when a defect is greater than .6 milli-inches (mils) deep. We want to be able 
to see defects this deep with great accuracy. 

The Michelson Interferometer consists of a light source that travels to 
a beam splitter and evenly splits the power of the source. Half the power of 
the source travels to a mirror that sometimes can be made to oscillate, called 
the reference arm, and the other half travels to the sample surface, called 
the sample arm. After the source is reflected off of the reference mirror and 
the sample, the light waves recombine at the beam splitter and are sent to a 
detector (see Fig. I). 1 
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Figure 1: Michelson Interferometer 
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When light travels through the Michelson Interferometer, the light goes 
to both the reference arm and the sample arm where each of the beams con- 
tains the same frequencies as the original light source. Assuming that the 
two paths are of equal distance and the two mirrors in the system reflect 
100% of the light, the light beams will recombine at the beam splitter and 
all the light from the source will go to the detector. When the mirror in the 
reference arm moves at a given frequency and amplitude through the point of 
equal distance, the two light paths will have total destructive and constructive 
interference with each other. This interference can be displayed on an oscillo- 
scope where the peak to peak voltage amplitude that is displayed represents 
the total alternating light intensity that is hitting the detector. When the os- 
cilloscope displays an oscillating voltage, rather than a constant voltage, we 
know that the waves from both paths are interfering and are thus coherent at 
that point. 2 When the oscilloscope is displaying a constant voltage, the light 
from both the paths are not interfering, meaning they are not coherent. To 
find coherence length, we can start where the amplitude of the light hitting 
the detector is at a maximum, and then move back the reference arm until 
the peak to peak amplitude has dropped by a half. Now, move the reference 
arm back through the maximum amplitude and stop where the amplitude 
has again dropped by a half. This total distance from half the maximum 
amplitude one way to half the maximum amplitude the other way (also called 
the full width at half maximum method, FWHM) is the coherence length. 2 
Since light waves can be added up linearly, the mathematics of the Michelson 
Interferometer will be the same if we send one frequency or many frequencies 
of light. Thus, if we pick a source that includes many frequencies of light, 
also known as a broadband light source, then the coherence length will be 
very small because each frequency in the source must be coherent and that 
only happens when the two paths of the Interferometer are lined up within a 
few microns. 3 

If we replace a sample with a defect in the sample arm, we can scan 
across the surface of the sample with our light source, and when the intensity 
of the light hitting the detector decreases, we know that there is a defect. 
Then we can focus into the defect by pulling the reference arm back which 
makes the light travel further into the defect. The problem is that in order to 
even see an amplitude change when scanning across a sample, the coherence 
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length of the source must be smaller than the defect depth. Otherwise, the 
source will be coherent with the surface of the sample, and when it passes 
over the defect, it will remain coherent with the bottom of the defect. Thus, 
the smaller the coherence length, the more accurately we can measure the 
depth of a defect. Because white light is very broadband and has a small 
coherence length, we will use it as our source. 

Mathematics of the Interferometer 

The wave equation denotes the propagation of light and is given by: 4 

d 2 u _ 2 d 2 u 
dt 2 ° dx 2 

Where 

c = speed of light 

Let u(t ) be a function that describes the amplitude, or electrical field 
of a given light source. The solution of the wave equation, u(t), can be made 
into a complex function, but for this project we will treat it as a real valued 
function. 

Power, p(t), is related to amplitude with the following proportionality: 

p(t) = pu 2 (t) 

Where 

(3 is the proportionality constant 

Now, when we are dealing with a Michelson Interferometer, the power 
of a light source is evenly distributed to each arm of the Interferometer. Thus, 
the electrical field down each arm, the reference arm u r (t ) and the sample arm 
u s (t), is related to the source’s electric field, u(t) by: 

1 

= —j=u 
\/2 

u s (t) = 

However, the light must travel from the arms back to the beam splitter, 
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so the amount of amplitude that makes it to the detector, Ud(t), if both arms 
are the same distance from the beam splitter, is: 

Ud{t ) = ^ e r u(t ) + ^e s u(t) 

Where 

el = fractional amount of energy transmitted through the reference arm 

el = fractional amount of energy transmitted through the sample arm 

The detector cannot respond to the actual optical frequencies and in- 
stead takes a time average from an arbitrary starting position To to T+ To. In 
other words, T » Also, the arms of the Interferometer are not the same 
distance away from the beam splitter because the reference arm is translating 
back and forth through the point of equal distance, so there is a time delay, 
r, between the two waves. The detector takes the power, p(t), from the light 
source and converts it to a current. This current is then converted to voltage 
using a known resistance by Ohm’s Law. Thus, the detector is reading volt- 
age, V'(To), which can be modeled by the following equation, assuming that 
the sample has a single reflection: 

1 r i r X 1 

V(T) = — Jt ° 0(-e r u(t - r r - At sin (ut + (j > )) + -e s u(t - T s )fdt 

Where 


t = continuous time that the detector is integrating over 
Tq = time response of the detector 

T r — time delay from the source to the detector through the reference arm 
t s — time delay from the source to the detector through the sample arm 
At = magnitude of the mirror oscillations in the reference arm 
uj = angular frequency of the mirror in the reference arm 
<p = phase of the mirror in the reference arm 
T = perceived time 
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The magnitude is squared so that the integrand we are taking is power. 
But, what if the sample was replaced with a sample with multiple reflections? 
The new formula would be the following: 5 


V ( T ) = Jr Jr +T ° P\l € ru(t - T r - Ar sin (ut + (f>)) + E ^ Sn u{t - r s J] 

-tQ JI * n = 1 * 


*L l 




Where 


N — number of discrete reflections in the sample 
Now, expand the integrand and factor out the (3 to get: 

V (T) = jr It +T °^ 4 u2 ^ ~ Tr ~ Ar sin ( ut + <£))+ 


~ru{t — r r — At sin {uit + <f>)) E e sn w (*- r s„) + E E ^^u(£-r Si )u(t-r Si )]dt 


Since (3u 2 {t) = p(t) and Ar is very small, we can treat the first term 
of the integral as a constant: 

P = Y 0 r ~4 u2 ^ ~ T r~ AT Sin ( ut + 4>))dt 

Which gives us: 


v (T)=p+iir°j^ 


T 0 


N 

T r - At sin (ut + <f>)) Y e s „u(t - r s J+ 

n= 1 


E E - T Si )u{t - r s .)dt 

i=lj=l 4 

The last term of the integrand is an autocorrelation function. When the 
reflections from the sample are far apart, meaning they are several coherence 
lengths away, then the autocorrelation function is zero. Thus the only terms 
that survive from this double sum are when i — j: 

V(T) = P + S. f T+To ± u ( t - r r - Ar sin (ut + 4>)) E e s „u(t - r s J+ 

Jo JT 4 n=l 

E e fAt - r Si )dt 

i = 1 4 
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Here, the last term is a sum of different powers, which can be taken 
out of the integral and combined with the constant power P to make a new 
constant, P': 

V{T) = P' + ^ F +T ° y u(t -T r - At sin (ut + <£)) £ e Sn u(t - r Sn )dt 

The sum can be moved to the outside of the integral along with the 
constants. We can assume that the time interval that we are averaging over 
is very large since the time interval that the oscilloscope is detecting is much 
larger than the actual time that a single wave of light takes to cross the 
detector, so we can let T — ► oo. 

V (T) = P' + jE e Sn lim ^ / T+T ° u(t - r r - At sin (ut + <f>))u(t - T Sn )dt 
4 n= i To— >oo 1qJT 0 

Now, we can treat the resulting integral as the autocorrelation of an 
ergodic process, r(r Sn ). 

V(T) = P' + ^Y. « fc r(T fc - TV - At sin (ut + </>)) 

4 n =l 

We know from the Wiener-Khinchin theorem that if we take the Fourier 
Transform of the autocorrelation function, we get the power spectrum (the 
proof is in the appendix). Since the Fourier transform is a linear transform, 
we get the following: 

Let T n = t Sji — T r — At sin (ut + 0). 

T[V(T)\ = />' + ^ E e«.r(T„)]e-*-dT,, 

J—oo 4^i 

- V 27 rP'S(u n ) + ^ E e s n p(^n) 

4 n— 1 

Where 

<5(^n) = the Dirac Delta function 

However, the power spectrum p(u n ) is centered about the center fre- 
quency of the light source. We must center the power spectrum over the 
origin so that we can single out the envelope of the autocorrelation. We can 
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also use the power spectrum to see if our results make sense with our light 
source. 

Let coo = the center frequency. 

R e n 

= V27T P'5{u n ) + ^ E e s n p{{v n - Wo) + W 0 ) 

4 n= 1 

The shift by + uiq makes the high frequency part of the autocorrelation 
function. Thus, we will take the Fourier transform again and use the shifting 
property of the transform to get: 4 

= rW^P'SM] + ^ p( w „ - «„)] 

4 71 = 1 

= -p' + ^ E 

4 n=l 

Where T e( t t>) is the envelope of the autocorrelation function. Substi- 
tute in what r n is and expand the exponent: 

= P' + ^E e Sn e-^ o(Tr - T ^)e iwoArsin M+^r E (r Sn - r r - At sin (cot + </>)) 

4 n =l 

Since the modulation At does not greatly affect the envelope of the 
autocorrelation function, we can treat it as zero which gives us: 

— iu 0 (r r - T Sn ) MJ 0 ATsin(w*+</>)p ( _ 

S n u Hf \ ' S n T 

Because we are using real valued functions, we will take the real part 
of this function to get: 

(5e r N 

= P' + “T E e s n cos (woC'Tr - T Sn )) cos (woArsin (cut + 0))r£;(r Sn - T r ) 

4 T>.= 1 


= p' + 


/?e, 


N 

E € 

77 — 1 


Where 

cos (wo(T r — r Sn )) = the high frequency term of the autocorrelation function 
cos (woArsin (2irft + (b)) = the modulation of the mirror 
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The setup that we have here at Kennedy Space Center involves a oscil- 
loscope that displays the voltage of the light source that is hitting the detector 
versus time. This oscilloscope is then connected to a spectrum analyzer which 
takes the Fourier transform of the voltage versus time graph and shows the 
voltage versus frequency graph. This enables us to look at the harmonics of 
the light source of a sample. The spectrum analyzer is then connected to a 
lock-in-amplifier which allows us to take the voltage off of a single harmonic, 
rather than averaging over a series of harmonics. 

Results 

We started with a superluminiscent diode (SLD) with a center fre- 
quency of 799 nm and a spectral width of 41 nm. We calculated the coher- 
ence length to be close to 11 microns. We placed a sample with a known 
defect into the sample arm and focused into the defect by pulling back the 
reference arm. Then we recorded the oscillating portion of the voltage on the 
oscilloscope in 10 micron increments, and obtained the graph in Fig. 2. 
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Figure 2: Defect Measurement with a SLD 

The high voltage at 9.8 mm shows the reflection off the front surface, 
but as the reference mirror is pushed back, we start to see the smaller inner 
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reflections from the defect at about 9.9 mm. From this data, we can deduce 
that the defect is around .1 mm deep. 

Now that we know we can see a defect, let us try to find smaller de- 
fects. The best way to do this is to use a light source with a small coherence 
length. White light has a broad bandwidth, and thus has a small coherence 
length. 3 When we replace a SLD with a white light source, the Interferometer 
becomes much more fragile in that it is very hard to find interference. Once 
interference is found, the coherence length can be calculated, again we will 
use the FWHM method again. The data that we obtained from a white light 
source by reading an oscilloscope and moving the reference arm in half micron 
increments away from the beam splitter is in Fig. 3. 
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Figure 3: White Light Coherence 

After fitting a gaussian to this data, we can use the FWHM method 
to find the coherence length. This data shows that we have a source with 
a coherence length of about 2.3 microns which is about .09 mils. Having a 
coherence length that is much shorter than our .6 mil cutoff point makes it 
easier to distinguish between the front reflection off the main surface of the 
sample and the defect reflections. 

With all the mathematics of the system laid out, let’s try to find the 
width of a thin glass cover slip by finding the point where the reference arm 
is coherent with the front surface, moving the reference arm back to find the 
coherent reflection off the back surface. We first measured the glass cover to 
be about 7 mils thick. This is still a factor of ten greater than the minimum 
depth a problem defect can be, but it will show us how easily this back reflec- 
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tion can be seen on the graph. Next, we replaced the 100% reflection mirror 
with the 4% reflection glass cover slip in the sample arm. Fig. 4 shows the 
data. 
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Figure 4: White Light Coherence with Two Reflection Sites 

The peak at around 8.96 mm shows the reflection off the front surface. 
The peak at around 9.23 mm is the reflection off the back surface. This shows 
that the front and back reflections are about 268 microns apart. When we 
convert 7 mils to microns, we get 177 microns. We also have to take into 
account the index of refraction in glass, which is about 1.5. Multiplying this 
index of refraction by the width of the glass, we get 265.5 microns, which is 
extremely close to our measured width of 268 microns, but that was not exact. 
The amplitude of the back reflection should only be 96% of the amplitude of 
the front reflection because the glass cover slip only has a 4% reflection. The 
probable reason for why the back reflection’s amplitude is much smaller than 
the front reflection is that the two surfaces were not exactly parallel because 
cover slips are not manufactured for this kind of experiment. 

Conclusion 

The short coherence length of white light enables us to distinguish 
very easily from reflections from the front and back of a thin glass cover slip. 
Looking at the last graph in the results section, the coherence from the front 
reflection falls very fast and allows for a wide gap between the end of the first 
reflection and the beginning of the back reflection. It would be very easy to 
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see a peak in the graph that is 15 microns from the peak of the front surface, 
which would enable us to see defects to the depth of the orbiter window’s min- 
imum problematic defect depth with great accuracy. To calculate the depth 
of a defect, find the front reflection surface and then move the reference arm 
back slowly to find the reflection from the bottom of the defect and see how 
far the reference arm moved back. Thus, this approach allows for an accurate 
and time-efficient way of measuring small defects in the orbiter windows as 
opposed to just detecting defects. 

Other Projects with the Interferometer 

Implementing the linearity of the Michelson Interferometer, we wanted 
to see if the coherence length would shorten when we added sources to the 
Interferometer. We thought that since white light is broadband, we could 
add SLD sources to get a broadband source that could theoretically be better 
than white light to find small defects. The next step was to see if the math- 
ematics worked in our favor. Assuming that the power spectrum of a light 
source is a gaussian or a sum of gaussians, we attempted to add several known 
power spectrums together. We noticed that if we added different spectrum 
together that overlapped each other, this would add to the spectral width of 
the spectrum and thus make a shorter coherence length. For example, we 
took a SLD with a center wavelength of 799 nm and a spectral width of 41 
nm and calculated a coherence length of about 13.4 microns. We then added 
another SLD with a center wavelength of 837 nm and a spectral width of 41 
and a coherence length of 23.41 microns. Fig. 5 shows the addition of these 
two spectrum which generates a wider spectral width. 

After shifting the gaussian to the origin and then taking the Fourier 
transform, we get the coherence function in Fig. 6. 

Using the FWHM method, we get a coherence length of 9.5 microns 
which is shorter than the shortest coherence length source. However, this still 
is much greater than the white light coherence length and it appears that we 
would have to add several sources to get the coherence length down to the 
coherence length of a white light source. Theoretically, we could add enough 
sources to the Michelson Interferometer to get a shorter coherence length 
than white light, but there are further problems. First, the power spectrums 
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Figure 5: Added SLD’s Frequency Spectrum 
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Figure 6: Coherence Function of Added Sources 


of each source must overlap the other power spectrums. If the power spec- 
trums do not overlap, then the sum of them will not affect the spectral width 
of either spectrum. Thus, we must have SLD’s with enough difference in 
central wavelength and spectral width to still overlap, but not lie right over 
each other (If they had similar central wavelength, then the spectral width 
would be the average of the two). Another problem involves the number of 
beam splitters required for a summation of light sources and the amount of 
power lost in each beam splitter. In order to add sources, Fig. 7 shows the 
modification to the Michelson Interferometer that could be applied. 

A new beam splitter would have to be added for every additional source. 
This would cause half the power of each source to be wasted. Thus, if we add 
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Figure 7: Michelson Interferometer with Two Sources 


5 sources we would lose 2 5 of the power from the furthest source from the 
detector. There are ways to prevent this huge loss in power, but the num- 
ber of beam splitters are still a problem. Ultimately, we want to be able to 
manufacture a hand-held device, and more beam splitters will only make it 
bigger. We concluded from this that adding sources would not be as efficient 
as simply using white light interferometry. 

The last thing that we tried, but did not finish entirely, was chang- 
ing the detector of the Interferometer to a spectrometer. Spectrometers have 
thousands of different sensors where each sensor corresponds to a small range 
of wavelengths and the spectrometer displays the intensity from a source hit- 
ting it. The benefit of this detector is that it does not require the reference 
arm of the Interferometer to oscillate. Therefore, there are no moving parts 
in this setup, which is a big plus. When we block one of the arms of the 
spectrometer with white light as our source, we obtain the spectrum in Fig. 
8. 

Also, when the two arms of the Interferometer are equal, we obtain a 
similar spectrum with a little more power. When we change the reference 
arm’s distance, we get oscillation in the spectrum. Assume that a single 
wavelength of light travels through the Interferometer, which can be modeled 
by: 

A cos ( kz — ut ) 


Where 


A = amplitude of the wave 

k — — and A = wavelength 
A 
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2 = distance the light wave traveled 
v = angular frequency of the light wave 
t = time 



Figure 8: Spectrum of White Light 

When we place this wavelength in the Michelson Interferometer, we get 
two power terms which is the DC offset, but what we are concerned with is 
the cross term: 


A 2 cos (kz refere nce - vt) cos (kz sample - vt) 

Where the subscripts on the z’s denote which path they travel down 
(refer back to the Mathematics of the Interferometer section). Using a trigono- 
metric identity, we can replace the cross term with: 

A 2 A 2 

— cos (deference — ^sample)) COS ^ ^sample ^reference) — 

The last term is a high frequency term, which we can ignore. This 
leaves us with a difference in distances from each arm: 

A 2 

— cos ( kAz )) 

Z 

Where 

Az = ^reference — ^sample 
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This shows that when the path lengths in the Interferometer are differ- 
ent, the number of peaks shown on the spectrum is about double the number 
of microns that the path lengths are actually different. It is double because 
the round trip of the light is taken into account. However, we can get more 
accurate measurements than just counting peaks because we can count frac- 
tions of peaks, thus we can see fractions of microns. 

We moved the reference arm back 5 microns from where the two paths 
of the Interferometer were equal and we obtained the spectrum in Fig. 9. 



Figure 9: Spectrum of White Light with different path lengths in the Interferometer 

Notice that the basic curve of the graph is still the power spectrum of 
the source. By dividing out the original spectrum of the source, we are left 
with the oscillation. What this comes down to is that when a wave is split in 
an Interferometer, the phase of each split wave is related to the difference in 
distance each wave travels. 

If we had more time we would be able to pursue this project further, 
but unfortunately, we ran out of time. Hopefully this idea can be devel- 
oped further and actually tested to see if it can actually find small defects in 
windows. 
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Appendix 


Proof of the Wiener-Khinchin Theorem 

The autocorrelation function r(r n ) is defined as 

/ OO 

u(t)u{T n + t)dt 

-OO 

The Fourier transform of u(t) is 

u(t) = / u{uj)e~ luit dw 

The Fourier transform of the complex conjugate is 

u(t) — [ u^e^duj 

Plug these into the autocorrelation function to get 

T£;( r n) = / [/ u{u)'e lu ' t dw']\l u(u>)e~ lu! ^ Tn+t ^duj]dt 

J —oo J —oo J— OO 

= f f f u(u')u((v)e~ ^t ^~^^e~ ^UTn dtdlJJduJ , 

J — oo J— oo J — oo 

= [ f u(u')u(uj)6(uj' — (jj)e~ lUTn du>dw' 

J — oo J —oo 

= j u(u)u(<jj)e~ tu>Tn duj 

J— OO 

— f \u(u))\ 2 e~ lUTn du 

J —OO 

= F[\uH\ 2 ](T n ). 

Therefore, the autocorrelation function is the Fourier transform of the power 
spectrum . 2 
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